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ABSTRACT 

We present VLT-FORS1 spectra of a sample of 34 faint 20.0 < g* < 21.1 A-type 
stars selected from the Sloan Digital Sky Survey Early Data Release, with the goal of 
measuring the velocity dispersion of blue horizontal branch (BHB) stars in the remote 
Galactic halo, R ~ 80 kpc. We show that colour selection with 1.08 < u* — g* < 1.40 
and —0.2 < g* — r* < —0.04 minimises contamination of the sample by less luminous 
blue stragglers. In classifying the stars we confine our attention to the 20 stars with 
spectra of signal-to-noise ratio > 15 A -1 . Classification produces a sample of eight BHB 
stars at distances 65 — 102 kpc from the Sun (mean 80 kpc), which represents the most 
distant sample of Galactic stars with measured radial velocities. The dispersion of the 
measured radial component of the velocity with respect to the centre of the Galaxy 
is 58 ± 15km s . This value is anomalously low in comparison with measured values 
for stars at smaller distances, as well as for satellites at similar distances. Seeking 
an explanation for the low measured velocity dispersion, further analysis reveals that 
six of the eight remote BHB stars are plausibly associated with a single orbit. Three 
previously known outer halo carbon stars also appear to belong to this stream. The 
velocity dispersion of all nine stars relative to the orbit is only 15 ± 4km s _1 . Further 
observations along the orbit are required to trace the full extent of this structure on 
the sky. 
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1 INTRODUCTION 

The existence of a dark massive halo appears to be a generic 
feature of many galaxies, but the total masses, sizes, and 
the formation history of galactic halos are poorly under- 
stood. This is mostly because we do not have large enough 
samples of dynamical tracers at sufficiently large radii. The 
formation and extent of such mass distributions are of great 
importance in understanding the nature of the dark matter 
and its role in galaxy formation and evolution. For instance, 
the quantification of the dark matter content of the Galaxy 
would allow us to construct a picture of the assembly of 
the various baryonic components, through comparison with 
simulations. The baryonic components can, in turn, provide 
information about the evolution of the halo. 
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There is compelling evidence that at least part of the 
stellar halo has been built up via the accretion of smaller 
satellite galaxies. Numerous searches have been made for 
streams of material responsible for building up the Galaxy. 
A striking example is the identification of the Sagittarius 
dwarf galaxy (Ibata, Gilmore & Irwin, 1994) and its stellar 
stream (e.g. Helmi & White, 2001). Recently, an extensive 
stream of stars has been uncovered in the halo of the An- 
dromeda galaxy (M31), revealing that it too is cannibalis- 
ing a small companion (e.g. Lewis et al. 2004). Such streams 
yield crucial information on the accretion history and forma- 
tion of galaxy halos. Extended stellar streams have also been 
used to constrain the mass of the Galactic halo (e.g. John- 
ston et al. 1999) and in M31 (Ibata et al. 2004). A number 
of authors have noted the possible evidence for streams in 
the distribution of the intrinsically rare, but very luminous, 
carbon stars. Sanduleak (1980) proposed the association of 
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a single carbon star with the Magellanic Stream and Tot- 
ten & Irwin (2000) made the general observation that the 
non-uniform distribution of carbon-stars in their extensive 
survey of the halo may indicate that a number of the stars 
are associated with streams. 

We have previously argued (Clewley et al., 2002, here- 
after Paper I) that blue horizontal branch (BHB) stars are 
an ideal population for exploring the outer reaches of the 
Halo. Like carbon stars they are luminous standard candles 
but are also far more numerous. 

BHB stars are A-type giants. A-type stars in the Galac- 
tic halo are easily identified, as they lie blueward of the 
main-sequence turnoff (e.g. Yanny et al., 2000, hereafter 
Y2000). Unfortunately assembling clean samples of remote 
r > 60kpc 1 BHB stars is made difficult by the existence 
of a contaminating population of high-surface-gravity A- 
type stars, the blue stragglers, that are between one and 
three mag. fainter. Previous analyses required high signal- 
to-noise ratio (S/N) spectroscopy to reliably separate these 
populations (e.g. Kinman, Suntzeff, and Kraft, 1994), mak- 
ing identification of BHB stars in the distant halo unfeasible. 

This paper is the third in a series. In Paper I we de- 
veloped two classification methods that enabled us to over- 
come the difficulties in cleanly separating BHB stars from 
blue stragglers, and outlined an observational programme to 
survey the halo for BHB stars. In the second paper (Clewley 
et al., 2004, hereafter Paper II), we presented photometry 
and spectroscopy of faint 16.0 < B < 19.5 candidate BHB 
stars in two northern high Galactic latitude fields and four 
southern fields. This work resulted in a sample of 60 BHB 
stars at distances 11 < R < 52 kpc (mean 28 kpc), with mea- 
sured radial velocities. Here we apply the methods of Papers 

1 and II to survey for halo BHB stars at much greater dis- 
tances, 65 < R < 115 kpc. The new survey uses Sloan Digital 
Sky Survey (SDSS) photometry to isolate a sample of faint 
halo A-type stars. Reliable classifications are derived from 
medium resolution spectroscopy using FORS1 at the VLT. 

The candidate BHB stars observed at the VLT were se- 
lected using u* g*r* photometry from the SDSS Early Data 
Release (EDR) data set (Stoughton et al., 2002). The EDR 
photometry was preliminary, and has since been revised. Be- 
cause we need accurate Johnson-Kron-Cousins B, V magni- 
tudes for the classification, we have used the more recent 
SDSS Data Release 2 (DR2) photometry (Abazajian, 2004) 
of the EDR-selected candidates for this purpose. In Section 

2 of this paper we describe the selection of the BHB candi- 
dates from the EDR data set, and provide our prescription 
for transforming the DR2 g,r magnitudes of A-type stars 
to B,V magnitudes. Section 3 provides a summary of the 
VLT spectroscopic observations, and the data reduction and 
line measurement procedures followed. In Section 4 we use 
the methods of Papers I and II to classify these stars into 
categories BHB and blue straggler, and provide a summary 
table of distances and radial velocities of the eight stars clas- 
sified as BHB. We compute the velocity dispersions of the 
two populations and compare them with previous work. In 
Section 5 we discuss the kinematics of the BHB stars. We 



1 In this paper we use the coordinate r to denote Galactocentric 
distances and the coordinate R to denote heliocentric distances. 
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Figure 1. Upper: Two colour plot of stellar objects in SDSS 
stripe 10, in the RA range 200° < a < 230°, the region includ- 
ing the Sagittarius tidal stream. Filled circles mark objects in 
the apparent magnitude range 18.8 < g* < 19.2, expected to 
be predominantly BHB stars. Open circles mark objects in the 
apparent magnitude range 20.5 < g* < 21.5, expected to be pre- 
dominantly blue stragglers. The box marks our adopted colour 
selection. The dashed line marks the division between high and 
low gravity A stars employed by Y2000. Lower: Histograms of the 
u* — g* colours of the stars in the upper plot. The thick line cor- 
responds to the filled circles, and the thin line to the open circles. 
Typical colour errors for these two samples are cr(u* — g*) = 0.05 
and 0.20, respectively. 

perform an orbital analysis of the sample and suggest that 
most of them are plausibly associated with a single orbit. 



2 SELECTING THE BHB CANDIDATES 
2.1 Colour selection 

We selected candidate BHB stars using the SDSS EDR point 
spread function (PSF) u*g*r* magnitudes of stellar objects. 
The SDSS photometry has evolved between EDR and DR2 
for a variety of reasons: i) the reference photometric system 
is now that of the SDSS 2.5m telescope itself, rather than 
the photometric monitoring telescope, ii) the absolute cali- 
bration of the standards has improved, and the calibration 
of the survey data relative to the standards has improved, 
iii) instrumental systematics (especially scattered light) are 
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now better understood. Therefore, we have taken advan- 
tage of the improved photometry in DR2 for the subsequent 
analysis, in particular in classifying the objects and esti- 
mating their distances. To distinguish between photometric 
systems, the EDR system is designated by asterisks, the sys- 
tem of the photometric monitoring telescope is designated 
by primes, and the DR2 magnitudes are unadorned, which 
is the SDSS convention. Bearing in mind that all the stars 
are in the remote halo, all the SDSS magnitudes discussed 
in this paper have been corrected for Galactic extinction, 
using the map of Schlegel, Finkbeiner & Davis (1998). 

We limited ourselves to the northern equatorial stripe 
in the EDR data set, which covers 145° < a < 236°, 
-1.25° < S < +1.25° (J2000). This is SDSS stripe 10, 
observed in runs 752 and 756 (Stoughton et al., 2002). In 
selecting candidate distant BHB stars from the EDR, we 
were guided by the results of Y2000, who studied the spa- 
tial distribution of a sample of A-type stars selected from 
the EDR using the (reddening-corrected) colour selection 
box -0.3 < g* - r* < 0.0, 0.8 < u* - g* < 1.5. In plot- 
ting apparent magnitude against a for these A-type stars, 
Y2000 discovered that the distribution is not smooth, but in- 
cludes striking over-dense regions occurring in bands. These 
have subsequently been identified as tidal debris from the 
Sagittarius dwarf galaxy. Furthermore the bands occur in 
pairs, coincident on the sky, but separated by ~ 2mag. The 
brighter objects are the more luminous BHB stars, and the 
fainter objects are blue stragglers at the same spatial loca- 
tion in the halo. 

Our aim was to use the EDR data set to select a sam- 
ple of candidate faint BHB stars, with minimal contamina- 
tion by blue stragglers, in order to make the most efficient 
use of the spectroscopic time awarded (for classification and 
velocity measurement). To reach large distances, we chose 
the magnitude range 20.0 < g* < 21.1, corresponding to 
~ 65 < R <~ 115 kpc, if the objects are BHB stars. For- 
tunately, as demonstrated by Lenz et al. (1998) using syn- 
thetic photometry, the u*g*r* colours of A-type stars show 
some dependence on surface gravity. Given the accuracy of 
the photometry at the distances of interest, ~ O.lmag., the 
u*g*r* colours cannot provide reliable separation of the two 
populations. Nevertheless Y2000 demonstrated that a colour 
cut in the u* —g* verus g* — r* plane is effective in enhancing 
the contrast of the individual bands of tidal debris i.e. can 
substantially reduce the contamination of one population by 
the other. 

Another way of looking at this is illustrated in Fig. 
1. The upper diagram plots the colours of all EDR stars 
within the colour selection box of Y2000, for the limited 
range 200° < a < 230°, which is the region of the EDR 
containing the strongest tidal debris bands. Stars in the 
brighter band 18.8 < g* < 19.2, which should be predom- 
inantly BHB stars, are marked with solid symbols, while 
stars in the fainter band 20.5 < g* < 21.5, which should 
be predominantly blue stragglers, are marked with open 
symbols 2 . For reference the dashed line shows the dividing 
line used by Y2000. The brighter band, typical colour er- 
ror a(u* — g*) — 0.05, is mostly confined to a narrow range 



2 The larger magnitude interval selected is because blue strag- 
glers have a larger spread in luminosity than BHB stars (§3.2). 



in u* — g* , which evidently defines the colour domain of 
the BHB stars. The fainter stars, open symbols, are con- 
centrated towards the top of the plot, but are spread over 
a larger colour range. However, much of the spread is ac- 
counted for by the larger colour errors a(u* — g*) = 0.2, as 
can be seen by reference to the histogram in the lower plot. 
It is evident that the mean u* — g* colour of blue stragglers 
is substantially bluer than for BHB stars. 

The spectroscopic classification criteria (detailed in §4) 
work best near (B — V)o = 0.1, which corresponds to 
g* — r* = —0.125, using the colour transformation provided 
by Fukugita et al. (1996). On this basis we adopted the 
colour cuts shown by the box in Fig. 1, defined by 1.08 < 
u* — g* < 1.40, —0.2 < g* — r* < —0.04, and limited candi- 
date selection to objects with colour error a(g* —r*) < 0.07. 
While these colour cuts should be nearly optimal in terms 
of the fraction of candidates that are BHB stars, we would 
still expect substantial contamination by blue stragglers, on 
account of the large u* — g* colour errors at the faint mag- 
nitudes of the sample, 20.0 < g* < 21.1. 

The distribution in a and g* of all the stars satisfying 
these criteria is shown in Fig. 2. The higher density of points 
at a > 200° is due to blue stragglers in the Sagittarius tidal 
stream. Our classification methods produce samples of BHB 
stars that are contaminated by blue stragglers at the level 
of < 10%, for samples of A-type stars with a typical mix 
of the two populations (Paper I). The contamination of our 
BHB sample would be substantially greater if we attempted 
to classify stars in this region, so we confined our sample to 
a < 200°. The clump visible at a = 153°, g* ~ 20.2, is the 
horizontal branch of the Sextans dwarf spheroidal 3 , centre 
a = 153.3°, 8 — —1.61°, J2000 (Irwin and Hatzidimitriou, 
1995). In order to avoid stars in Sextans we confined our se- 
lection to a > 160°. The final selection includes 54 objects, 
of which we observed the 35 listed in Table 1. One of these, 
no. 32, proved to be a quasar. Looking ahead, of the remain- 
ing 34 candidates we are able to classify 20. These classifica- 
tions are indicated on Fig. 2, with solid circles representing 
the eight stars classified BHB, and open circles representing 
the 12 stars classified blue straggler. For the remaining 14 
candidates the classifications are uncertain, mostly due to 
insufficient S/N. 

In Table 1 we provide details of the 35 objects observed. 
Column 1 is our running number, and column 2 lists the co- 
ordinates. Successive columns provide the dereddened DR2 
g magnitude, and the dereddened u — g and g — r colours. 
The last column provides the dereddened B — V colour, cal- 
culated using the transformation derived in §2.2. There is 
good agreement between the EDR and DR2 photometry in 
the mean, but with noticeable scatter. For example, looking 
at the difference (g* —r*) — (g — r) for our targets, the mean 
is 0.00 mag., and the standard deviation is 0.03 mag. 

2.2 Transforming from g — r to B — V 

We use two methods to classify stars into categories BHB 
star and blue straggler. One method makes use of (B — 
V)o colours. Therefore we need to convert the extinction 

3 Note that the discussion in Y2000 of these objects mistakenly 
cites a paper concerned with a different galaxy, Sextans A. 



4 L. Clewley et al. 



No. 


Identification (J2000) 






n 

y 








(u 




a) 
J) 





(a - 


1 n> 


(B 




V 


In 


01 


71 0531 9 


.31 




.79 


20. 


.09 


± 





.02 


1 


.08 


± 





09 


-0.08 


± 0.03 


0. 


15 


± 





,03 


02 


T1 05322 


38 


— 004449 


.73 


20 


.19 


± 





.03 





.81 


± 





,16 


-0.10 


± 0.05 


0. 


14 


± 





.01 


03 


T1 1 0639 


10 


— 004720 


06 


20 


16 


± 





.03 


1 


11 


± 


0. 


,15 


-0.13 


± 0.05 


0. 


11 


± 


0. 


04 


04 


T1 1 93^9 


91 


— 00^71 Q 


83 


20 


01 


± 





.02 


1 


33 


± 


0. 


11 


-0.16 


± 0.03 


0. 


09 


± 


0. 


03 


05 


J112744, 


35 


_i_ 0D1 "SOS 


.47 


20 


.85 


± 





.03 


1 


.29 


± 





,18 


-0.14 


± 0.05 


0. 


11 


± 





.01 


06 


J -L lOUOU . 


42 


— 005147 


\\ 


20. 


.27 


± 





.04 


1 


.03 


± 





,10 


-0.14 


± 0.05 


o 


\\ 


± 





,04 


07 


T1 1 ^zL07 


58 


— 0047^ 


65 


20 


21 


± 





.04 


1 


.28 


± 


0. 


19 


-0.06 


± 0.06 


o 


18 


± 


0. 


04 


08 


T1 1 481 8 


87 


— 00^Q91 

WOW J- 


.03 


20. 


.06 


± 





.05 


1 


.10 


± 





,15 


-0.08 


± 0.07 


0. 


16 


± 





,05 


09 


T1 1 5525 


99 


— 003601 

\J KJ <J KJ \ ) J. 


.66 


20. 


96 


± 





.03 


1 


.27 


± 





,24 


-0.08 


± 0.06 


0. 


16 


± 





,04 


10 


J120024. 


18 


4- 01 1026 


81 


21 


06 


± 





.04 


1 


.28 


± 


0. 


,24 


-0.13 


± 0.07 


0. 


12 


± 


0. 


05 


11 


T1 20855 


70 


4- 01 0929 


.20 


20 


.35 


± 





.03 


1 


.18 


± 


0. 


.12 


-0.14 


± 0.04 


0. 


10 


± 


0. 


,04 


12 


J121447. 


.45 


4- 00 100 1 


.27 


21. 


00 


± 





.04 


1 


15 


± 


0. 


17 


-0.10 


± 0.06 


0. 


14 


± 


0. 


05 


13 


T1 22305 


31 


— 011443 


.14 


20. 


.75 


± 





.04 


1 


.22 


± 


0. 


,21 


-0.10 


± 0.07 


0. 


14 


± 


0. 


,06 


14 


T1 9979R 


77 




03 


21 


.01 


± 





.03 


1 


.16 


± 





,18 


-0.09 


± 0.06 


o 


15 


± 





.05 


15 


T1 22802 


26 


— 01 Ql^l 


.91 


20 


93 


± 





.05 


1 


10 


± 


0. 


.22 


-0.11 


± 0.08 


0. 


13 


± 


0. 


06 


16 


T1 9^80^ 


94 




07 


21 


.04 


± 





.04 


1 


13 


± 


0. 


,22 


-0.06 


± 0.07 


o 


17 


± 


0. 


.05 


17 


T1 941 1 9 


31 


— 010447 


.12 


21 


.03 


± 





.04 


1 


.06 


± 





,22 


-0.06 


± 0.07 


0. 


17 


± 





,06 


18 


T1 24620 


96 


_ 009809 


.04 


20 


13 


± 





.02 


1 


26 


± 


0. 


,12 


-0.20 


± 0.04 


0. 


06 


± 


0. 


03 


19 


T1 24851 


75 


4- 003045 


.08 


20 


.05 


± 





.02 


1 


21 


± 


0. 


.10 


-0.07 


± 0.04 


0. 


17 


± 


0. 


03 


20 


T195004 


94 


4- 003422 


85 


20 


91 


± 





.04 


1 


.14 


± 





,22 


-0.10 


± 0.06 


0. 


14 


± 





05 


21 


71 25054 


28 


— 000759 

KJ KJ KJ 1 UiS 


.84 


20 


.44 


± 





.03 


1 


.12 


± 





,11 


-0.17 


± 0.04 


0. 


08 


± 





.04 


22 


71 251 23 


59 


4- 000345 


.91 


20 


.78 


± 





.03 


1 


.01 


± 





,15 


-0.11 


± 0.05 


0. 


13 


± 





.01 


23 


71 25308 


36 


— 000554 


.68 


20 


.75 


± 





.03 


1 


20 


± 


0. 


,17 


-0.08 


± 0.05 


0. 


15 


± 


0. 


04 


24 


J125336 


.47 


- 002414 


.48 


20 


.01 


± 





.03 


1 


.16 


± 


0. 


09 


-0.14 


± 0.04 


0. 


11 


± 


0. 


,03 


25 


J125934. 


17 


+ 002058 


.64 


21 


.07 


± 





.04 


1 


.00 


± 





,22 


-0.07 


± 0.07 


0. 


17 


± 





,06 


26 


J125938. 


97 


+ 000748 


.09 


21. 


.07 


± 





.04 


1 


.56 


± 


0. 


.32 


-0.10 


± 0.07 


0. 


14 


± 





,05 


27 


J130528. 


14 


+ 004855 


.90 


20. 


23 


± 





.04 


1 


,14 


± 


0. 


11 


-0.06 


± 0.05 


0. 


17 


± 


0. 


,04 


28 


J130907. 


25 


+ 005731 


.22 


20 


11 


± 





.03 


1 


16 


± 


0. 


11 


-0.10 


± 0.05 


0. 


14 


± 


0. 


04 


29 


J131109. 


31 


+ 000950 


.36 


20. 


.67 


± 





.03 


1 


.28 


± 


0. 


,20 


-0.18 


± 0.05 


0. 


07 


± 





,04 


30 


J131155. 


.98 


+ 003914 


.93 


20. 


61 


± 





.03 





.84 


± 


0. 


13 


-0.11 


± 0.05 


0. 


13 


± 


0. 


,04 


31 


J131252. 


50 


+ 000821 


.47 


20 


.72 


± 





.03 


1 


05 


± 


0. 


16 


-0.13 


± 0.05 


0. 


11 


± 





04 


32 


J131435. 


78 


+ 010329 


82Q 


20 


.55 


± 





.03 


0. 


.98 


± 


0. 


.16 


0.01 ± 0.06 


0. 


22 


± 


0. 


04 


33 


J131458. 


12 


- 004706 


.39 


20 


.30 


± 





.03 


1 


.11 


± 





,14 


-0.08 ± 0.04 


0. 


15 


± 





,04 


34 


J131538. 


01 


- 010853 


.30 


20 


.26 


± 





.03 


1 


.25 


± 





.13 


-0.08 ± 0.04 


0. 


16 


± 


0. 


,04 


35 


J131602. 


90 


+ 010150 


.60 


20 


50 


± 





.03 


1 


,24 


± 


0. 


,17 


-0.10 ± 0.05 


0. 


.14 


± 


0. 


04 



Table 1. Photometric data for the 35 BHB candidates from the SDSS DR2. One candidate that was subsequently found to be a quasar 
has 'Q' appended to the SDSS identifier. 



corrected g — r colours to B — V . Fukugita et al (1996) 
quote the transformation 

B — V = (g — r' + 0.23)/1.05, (1) 

computed from synthetic photometry. Smith et al. (2002) 
quote the transformation 

B-V = (ff'-/ + 0.19)/0.98, (2) 

based on observations of standard stars. These two relations 
are plotted in Fig. 3, as the dashed and dotted lines re- 
spectively. The two linear relations differ somewhat even 
over the narrow range of colours of interest in this paper, 
—0.2 < g—r < —0.04, by some 0.03 mag. Whereas one would 
naturally prefer the empirical relation over the synthetic re- 
lation, there is some indication that the actual transforma- 
tion is non-linear in the region of the A stars, as the four 
stars measured by Smith et al. (2002) in the colour range of 
interest, lie on average 0.04 mag. above the linear relation 
(which is a fit over a wide colour range). Further evidence 
that the relation is non-linear comes from our own photom- 
etry, which is less precise, but has many more stars. Plotted 
in Fig. 1 are the colours of the 60 stars in the colour range 
-0.05 < (B — V) < 0.40, with (B - V) measured by our- 



selves (Paper II), which also have SDSS DR2 g — r colours. 
Only the errors on (B — V)o are plotted, as these dominate. 
Our data points are systematically high relative to the lin- 
ear relation of Smith et al. (2002). An additional source of 
uncertainty is the fact that the above linear relations were 
derived for the u' g'r'i'z' system of the photometric moni- 
toring telescope, slightly different from the ugriz system of 
the SDSS 2.5m telescope itself, used for DR2. 

Accurate (B — V)o photometry is required for the classi- 
fication of the stars using the D0.15 -Colour method (§4). A 
systematic error in the (B — V)o colour as large as 0.05 
mag. could result in many of the classifications being in 
error. Therefore we re-investigated the colour transforma- 
tion. We computed synthetic colours using the methods de- 
tailed in Hewett et al. (in prep.). We used model stars with 
[Fe/H]= —1 and log <jr=3.5 (an appropriate surface grav- 
ity midway between BHB stars and blue stragglers), from 
Kurucz (1993). We fit a cubic polynomial to the relation be- 
tween the (B — V)o and (g — r)o synthetic colours, for the 
colour range of Fig. 3. Finally, bearing in mind the uncer- 
tainty in the absolute calibration of the SDSS magnitudes 
onto the AB system (Fukugita et al., 1996), we allowed the 
zero point of the relation to be a free parameter, established 
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Figure 2. A plot of a against g* for all the BHB star candidates 
in EDR stripe 10, satisfying our selection criteria (detailed in the 
text). Selection was confined to the range 160° < a < 200°. There 
are 54 stars satisfying the selection criteria, of which the 35 listed 
in Table 1 were observed and 20 were classifiable. The eight stars 
classified as BHB stars are marked by filled circles, and the 12 
stars classified as blue stragglers are marked by open circles. 



by shifting the derived curve vertically to give the best fit 
to the data of Fig. 3. The curve provides a better fit than 
the two linear relations plotted. Furthermore the average 
offset of the four stars measured by Smith et al. (2002) in 
the colour range of interest reduces to 0.005 mag. We have 
therefore adopted this relation. The transformation is given 
by 

B — V = 0.764( 5 -r)-0.170(g-r) 2 +0.715(3-r) 3 +0.218, (3) 

and is plotted in Fig. 3 as the bold solid line. We stress 
that this relation is specifically for A-type stars, and is not 
expected to be reliable for other types of star. The adjust- 
ment to the zero point was very small, only 0.014 mag. The 
agreement is encouragingly good and gives considerable con- 
fidence in the many elements going into this comparison - 
the SDSS DR2 photometry, the absolute calibration of Vega, 
and of the SDSS standard stars, our own photometry, the 
measurement of the different passband response functions, 
and the synthetic stellar spectra. The computed (B — V)o 
colours of our targets are listed in the final column of Table 
□ 

In summary, we have presented evidence that, in the 
colour range of interest, A stars lie systematically off the 
simple linear colour transformation measured by Smith et 
al. (2002), by about 0.05 mag. and we have derived a cubic 
relation that provides an improved fit. 



3 SPECTROSCOPIC OBSERVATIONS, 
ANALYSIS, AND RESULTS 

3.1 Observations 

We used the VLT FORS1 instrument, in service mode, 
over two periods, from 2003/03/25 to 2003/04/09, and from 
2004/01/30 to 2004/03/20, to obtain medium resolution op- 
tical spectra of the 35 BHB candidates. The instrument is 
equipped with a 2048 2 Tek CCD, with a projected scale of 




Figure 3. Colour transformation between dereddened SDSS DR2 
(g — r )o colour and (B — V)o- The data plot the photometry of 
the 60 stars from Paper II with SDSS photometry. The dashed 
line is the transformation derived by Fukugita et al. (1996) from 
synthetic photometry, and the dotted line is the transformation 
measured by Smith et al. (2002). The solid line is the cubic rela- 
tion derived here from synthetic colours of Kurucz (1993) model 
stars. The vertical lines mark the (g* — r*)o colour selection limits. 



0.2"pixel _1 . We used the 600B grating, giving a dispersion 
of 1.2 A pixel -1 . With the 0.7" slit, the resolution achieved, 
measured from arc lines, was about 4A, which is sufficient for 
the line-fitting procedure. The spectral coverage was 3400- 
5700 A, which includes the relevant lines H<5, H7, and Ca 
II K A3933A. Three BHB radial velocity standard stars in 
the globular cluster M5 were observed twice each. Table [5] 
summarises relevant information on the standards. Columns 
(1) to (5) list the identification, RA and Dec, V magnitude, 
(B — V)o colour, and the heliocentric radial velocity, Vq. 
The information in successive columns (6) to (11) in Ta- 
ble 3 contains averages of HS and Hy line measurements. 
Columns (6) to (8) in Table [5] list, respectively, the param- 
eters Do. 15, b, and c (explained in §4). The errors on the 
parameters b and c are provided in columns (9) to (11) in 
the form of A and B, the semi-major and semi-minor axes 
of the error ellipse in the b — c plane, and 6 the orienta- 
tion of the semi-major axis, measured anti-clockwise from 
the 6-axis. Here the error corresponds to the 68% confidence 
interval for each axis in isolation (see Paper I for further 
details). 

The requested integration times of between 765 and 
1980 seconds, for stars in the magnitude range 20.0 < g* < 
21.1, were estimated using the FORS1 exposure-time calcu- 
lator, on the basis of the requested seeing, transparency, and 
lunar phase, in order to achieve the minimum continuum 
S/N ratio of 15 A required to classify the stars (Paper 
I). All targets were observed near culmination, with a mean 
airmass of 1.20 ± 0.1. In the event, observations of 12 of the 
34 targets failed to achieve the required S/N, and therefore 
these targets cannot be reliably classified. The failures were 
primarily in the cases where the seeing was poor. In retro- 
spect, for these service observations we should have included 
a safety margin in the requested integration times. All im- 
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age frames were automatically bias and flat -field corrected 
by the FORS pipeline, and we then followed standard pro- 
cedures for sky subtraction, spectral extraction, and wave- 
length calibration. An error spectrum, used for the line pro- 
file fits, was computed from Poisson considerations. 

For the wavelength calibration, HgCdHe arc observa- 
tions were made during the day, and were used to derive 
the dispersion solution. To account for any flexure of the 
instrument, the zero point of the dispersion solution was 
established for each spectrum using the [OI] night-sky line 
at 5577.34 A. We found a rms drift of the zero point of 
fSkms" 1 over the entire data set. 

3.2 Analysis 

As stated earlier, one of the candidates proved to be a 
quasar. In the remainder of the paper we ignore this ob- 
ject, and refer only to the 34 stars. The spectra were used 
to measure the shapes and widths of the H<5 and H7 lines 
(for classification), the EW of the Call K line (to measure 
the metallicity), and the radial velocity of each star (for 
future dynamical analysis). For these measurements we fol- 
lowed the procedures set out in Papers I and II exactly, and 
we refer the reader to those papers for full details. Below 
we provide a brief explanation of how the Balmer lines were 
measured, and how the Call K line EW is used to deter- 
mine the metallicity. We then summarise how the magni- 
tude, colour, and metallicity are combined with the classifi- 
cation to estimate a distance for each star. 

Balmer line profiles. After normalising each spectrum 
to the continuum, we fit a Sersic function, convolved with 
a Gaussian of FWHM the instrumental resolution, to the 
H<5 and H7 lines. Two parameters of the fit, the scale width 
b, and the shape index c, are recorded. One classification 
procedure, the Scale width-Shape method, plots these two 
quantities against each other. A third quantity Do. 15, which 
is the line width at a depth 15% below the continuum, is 
derived from b and c. The second classification method, the 
D0.15 — colour method, plots this quantity against (5 — V)o- 
Because Do. 15 is a function of b and c, the two classification 
methods are not completely independent. 

Metallicities from Call K lines. The Call K line is the 
strongest metal line present in the wavelength range cov- 
ered by the spectra, and the only useful line in moderate 
resolution blue spectra for measuring metallicity. Plotting 
Call K line EW against (B — V)o, the metallicity is deter- 
mined by interpolation between lines of constant metallicity 
on this diagram (see Fig. 4). The uncertainty is established 
from the uncertainties of the two quantities plotted, and 
an additional uncertainty of 0.3dex is added in quadrature. 
This is the systematic error, and was established by com- 
paring metallicities derived by this method using high S/N 
data of comparable resolution, with accurate metallicities 
determined from high-resolution spectra. No attempt has 
been made to remove the possible contribution of interstel- 
lar Call K absorption from the stellar K measurements. For 
a remote halo star the typical Call K EW is 0.1lA/sin6 
(Bowen, 1991), with a 95% range of (0.06 - 0.31)A/sin6. 
This range translates to 0.07 — 0.35A for our fields. Our 
34 stars have mean EW 1.7A and standard deviation 0.8A, 
with only two stars having EW below 0.7A. Therefore inter- 
stellar Call K absorption is insignificant for the majority of 



our targets, but could bias the measured metallicities high 
for the small fraction of stars with the weakest lines. 

Distances. The absolute magnitude of a BHB star 
Mv (BHB) depends on both metallicity and colour (i.e. tem- 
perature). In Paper II we derived a relation for the abso- 
lute magnitude of BHB stars by combining published re- 
lations for the dependence of Mv on metallicity and on 
(B — Vq) colour, with the measured absolute magnitude at 
fixed metallicity and colour, as follows. The slope of the re- 
lation between apparent magnitude and metallicity for RR 
Lyrae stars in the Large Magellanic Cloud was measured by 
Clementini (2003) from observations of some 100 stars. We 
combined this with the measurement of the absolute mag- 
nitude of RR Lyrae stars at fixed metallicity, determined by 
Gould & Popowski (1998) from Hipparcos statistical paral- 
laxes, to derive the linear relation for RR Lyrae stars: 

M V (RR) = 1.112 +0.214[Fe/H]. (4) 

We then adopted a cubic expression determined by Pre- 
ston et al. (1991), for the (B — V)o colour dependence of 
the difference in absolute magnitudes between BHB and RR 
Lyrae stars, to produce the final expression for the absolute 
magnitude of BHB stars: 

Mv(BHB) = 1.552 + 0.214[Fe/H] - 4.423(5 - V)o 

+17.74(5- V)l- 35.73(5 -V)%. (5) 

Distances and associated errors are then determined 
using the apparent magnitudes Vb, and the corresponding 
photometric and metallicity errors. To compute V, we used 
the relation V = g' - 0.53(g' - r') (Fukugita et al., 1996), 
here disregarding the subtle differences between the differ- 
ent SDSS magnitudes (g,g',g*, etc.). The result produces 
distance errors of 6 — 10% for our confirmed BHB stars. The 
exact form of Equation particularly the zero point, re- 
mains controversial. Currently there are at least ten meth- 
ods of determining the absolute magnitudes of RR Lyrae 
stars. We refer the interested reader to a recent review of 
this subject by Cacciari & Clementini (2003). 4 

The absolute magnitudes of blue stragglers have been 
less well studied. Since we will not use the blue stragglers 
in any dynamical analysis, their distances are less interest- 
ing. In Paper II we adopted the following relation derived 
by KSK from data for blue stragglers in globular clusters 
published by Sarajedini (1993) 

Mv(BS) = 1.32 + 4.05(5 - V) - 0.45[Fe/H]. (7) 
3.3 Results 

The results of these measurements for the 34 candidate BHB 
stars are provided in Table|^| Column (1) gives the number 
of the star, and columns (2) and (3) record the EW of the H7 
line, and the spectrum continuum S/N per A. Our classifi- 
cation methods were developed specifically for objects with 
strong Balmer lines, defined by EW H7 > 13A, and with 
continuum S/N > 15A" 1 . In all, only 20 of the 34 candi- 
dates meet both criteria, and are therefore classifiable. The 

4 They find, by averaging over all ten methods in their review, 
M V (RR) = 0.93 ± 0.12 + (0.23 ± 0.04) [Fe/H] (6) 
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(B-V) o Power, c(yS) 

Figure 5. Classification of the 20 survey stars, using the Do. 15 -Colour (a) and the Scale width-Shape (b) classification methods. The 
solid curves are the classification boundaries explained in the text. Filled circles are stars classified BHB in both plots, i.e. the stars below 
the classification boundary in each plot. Filled triangles are stars below the classification boundary in only one plot but are nonetheless 
classified BHB, whereas open triangles are classified A/BS. In plot (a) there are eight stars below the boundary and 12 above it. In plot 
(b) there are seven stars below the boundary and 13 above it. A total of eight stars are classified BHB. 



ID 


RA (J2000) Dec. 


V 


(B - V) 


[km s" 1 ] 


-Do.i5(7< 5 ) 

[A] 


b( 7 S) 

[A] 


c(7<5) 


A 


B 


6 


M5-I-53 


J151836.35 + 020744.6 


15.06 


0.06 


52.2 ± 1.4 


29.47 ± 0.26 
29.38 ± 0.28 


8.22 

8.14 


0.91 
0.90 


0.11 
0.11 


0.014 
0.014 


1.511 
1.511 


M5-II-78 


J151826.93 + 020717.8 


14.95 


0.12 


42.2 ± 1.1 


27.81 ± 0.28 
28.22 ± 0.28 


7.58 
7.65 


0.87 
0.87 


0.11 
0.12 


0.015 
0.014 


1.509 
1.511 


M5-IV-05 


J151835.34 + 020227.9 


15.15 


0.15 


56.9 ± 1.2 


29.51 ± 0.28 
29.66 ± 0.27 


8.07 
8.12 


0.88 
0.89 


0.11 
0.11 


0.014 
0.013 


1.512 
1.512 



Table 2. Spectroscopic measurements of three M5 globular cluster BHB stars. The names are from Arp (1955) and Arp (1962), the 
photometry is from Cudworth (1979), and the radial velocities are from Peterson (1983). 



majority fail because of inadequate S/N . The information in 
successive columns (4) to (9) in Table 3 contain averages of 
H8 and line measurements. The quantities in columns 
(4) to (9) in Table 3, i.e. Do. 15, b, and c, A, B and 9, are 
also provided for the radial velocity standards in columns 
(6) to (11) of Table 2. Comparing these quantities enables 
us to use the observations of the standard stars as a further 
check of our classification methods. 

In column (10) of Tablets listed the EW of the Call 
K line. This is plotted against (B — V)o in Fig.|I] The mea- 
sured metallicity for each star, and the error (including ran- 
dom and systematic contributions) is provided in column 
(11). The large errors are a consequence of the compara- 
tively large errors in the g* — r* colours. The mean mea- 
sured metallicity of the stars plotted is —1.4 with dispersion 
0.6, similar to the mean value measured for our sample of 
brighter A-type stars (Paper II). There are no significant 
outliers, but this is not a strong statement, given the large 
errors. 

The radial velocity, corrected to the heliocentric frame, 
is provided in column (12), and the estimated distance, 



based on the classification from the following section, is pro- 
vided in column (13). 



4 CLASSIFICATION AND VELOCITY 
DISPERSION 

4.1 Classification 

As noted above, of the 34 candidates, only 20 meet the re- 
quirements on spectroscopic S/N and H7 EW for reliable 
classification. In the following we restrict our discussion to 
the classification of these 20 objects. Of the other 14 can- 
didates, five stars have EW H7 < 13A, and are considered 
unclassifiable. For the remaining nine candidates the S/N 
of the spectra is too low for the classification to be reliable. 
We have nevertheless followed the classification procedures 
for these objects, but for clarity have omitted them from 
Figures 4 and 5. The final classifications are flagged as ques- 
tionable. We have followed the classification procedures of 
Paper II (which are slightly different from those of Paper I) 
exactly, with the exception that we weight the two classifi- 
cation methods unequally, as detailed below. 
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EW( T | S/N 

[A] [A]" 1 



00.15(715) 

[A] 



!>(-?«) 

[A] 



c( T <5) 



EW(CallK) 
[A] 



[Fo/H] 



ob./cla 



[kpo; 



(1) 


(2 J 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 




(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


01 


19.2 


10 


3 


35.78±1.69 


9.24 





77 


0.68 


0.068 


1 


521 


0.92±0 


19 


-2.09±0.41 


86.6±19.0 


73.4±5.6 


(0.39) A/BS? 


02 


19.8 


15 


2 


41.98± 1.23 


11.12 





84 


0.47 


0.039 


1 


532 


0.34±0 


17 


~2.87±0.37 


393.8±12.4 


26.6±7.9 


(0.00) A/BS 


03 


17.7 


15 


9 


38.50±1.15 


10.95 





95 


. 4 6 


0.051 


1 


5 2 3 


1.14±0 


17 


-1.64±0.41 


260.9±12.7 


36.3±5.1 


(0.09) A/BS 


04 


20.0 


19 


2 


41.94±0.97 


11.27 





86 


0.38 


0.031 


1 


5 2 9 


0.50±0 


16 


-2.32±0.38 


172.0±10.9 


30.2±5.6 


(0.00) A/BS 


05 


16.1 


20.4 


30.97±0.81 


8.67 





90 


0.33 


0.043 


1 


507 


0.93±0 


15 


-1.78±0.45 


172.8±16.4 


98.5±7.9 


(0.97) BHB 


06 


18.8 


15 


7 


39.53±1.12 


10.81 





87 


0.44 


0.043 


1 


528 


1.44±0 


17 


-1.37±0.43 


354.3±13.8 


39.2±5.5 


(0.04) A/BS 


07 


19.1 


14 




38.75±1. 13 


9.44 





71 


0.45 


0.033 


1 


527 


1.35±0 


17 


-1.83±0.49 


-35.5±09.8 


30.3±6.8 


(0.03) A/BS? 


08 


18.8 


17.9 


38.38±1.14 


9.73 





78 


0.45 


0.037 


1 


526 


0.79±0 


17 


-2.29±0.51 


56.9±16.5 


27.0±7.1 


(0.08) A/BS 


09 


11 .2 


21 


2 


26.40±0.79 


5.52 





61 


0.30 


0.022 


1 


511 


1.88±0 


15 


-1.35±0.48 


45.0±09.4 


49.8±7.6 


(0.33) A/BS 


10 


20.4 


15.0 


43.93±1.18 


12.60 





96 


0.48 


0.047 


1 


527 


0.72±0 


17 


-2.11±0.55 


64. 5± 15.4 


47.8±9.9 


(0.00) A/BS 


11 


19.0 


12 


6 


39.75±1.49 


9.68 





74 


0.59 


0.040 


1 


528 


1.59±0 


18 


-1.20±0.42 


133.1±47.4 


43.0±5.0 


(0.05) A/BS? 


12 


1 1 . 6 


17 


5 


24.59±0.89 


5.71 





70 


0.34 


0.034 


1 


507 


2.55±0 


15 


-0.81±0.54 


96. 7± 15.2 


98.1±6.6 


unclassifiablc 


13 


16.8 


12 


3 


35.49±1.38 


9.14 





77 


0.55 


0.050 


1 


520 


2.18±0 


18 


-1.08±0.57 


160.2±17.6 


49.0±6.9 


(0.28) A/BS? 


14 


13.7 


21 


4 


29.05±0.75 


7.01 





73 


0.40 


0.040 


1 


514 


2.19±0 


11 


-1.14±0.48 


-21.2±10.5 


102.3±5.1 


(0.75) BHB 


15 


16.8 


9.0 


32.74±1.87 


7.95 





74 


0.77 


0.059 


1 


518 


2.86±0.20 


-0.50±0.71 


19.4±24.5 


91.4±6.5 


(0.52) BHB? 


16 


17.4 


15 


8 


36.82±1.08 


9.26 





76 


0.42 


0.036 


1 


526 


0.80±0 


16 


-2.36±0.49 


2.7±10.7 


40.3±10.9 


(0.15) A/BS 


17 


9.9 


14 


7 


21.03±1.17 


4.18 





57 


0.40 


0.042 


1 


492 


2.98±0 


11 


-0.80±0.61 


97.4±16.9 


101.5±7.3 


unclassifiablc 


18 


19.2 


19 


9 


39.32±0.91 


10.49 





85 


0.36 


0.032 


1 


526 


1.94±0 


16 


-0.49±0.46 


58.3±10.0 


49.4±3.3 


(0.01) A/BS 


19 


15.0 


16 


8 


33.42±1.02 


8.24 





76 


0.41 


0.033 


1 


522 


2.06±0 


16 


-1.31±0.41 


38.3±12.2 


68.7±4.8 


(0.65) BHB 


20 


18.4 


13 


1 


35.98±1.41 


8.63 





72 


0.56 


0.041 


1 


525 


1.28±0 


18 


-1.68±0.47 


120.2±13.4 


46.2±8.5 


(0.24) A/BS? 


21 


19.2 


15 


7 


40. 56 ±1.21 


12.02 





99 


0.47 


0.059 


1 


529 


0.95±0 


17 


-1.63±0.38 


-95. 6 ±20.0 


42.1±5.9 


(0.01) A/BS 


22 


20.2 


13 


2 


38.64±1.57 


8.84 


0.67 


0.60 


0.037 


1 


530 


3.13±0 


19 


-0.38±0.55 


74.2±22.5 


58.5±5.1 


(0.06) A/BS? 


23 


16.1 


13 


9 


34.92±1.33 


7.96 





68 


0.54 


0.035 


1 


525 


1.73±0 


18 


-1.44±0.45 


-31.5±16.6 


44.2±7.2 


(0.29) A/BS? 


24 


13.7 


21 


9 


31.28±0.75 


8.27 


0.84 


0.29 


0.033 


1 


516 


1.01±0 


15 


-1.71±0.36 


18.0±09.7 


66.8±5.0 


(0.96) BHB 


25 


12.2 


16 


1 


24.90±1.08 


5.55 





66 


0.41 


0.037 


1 


508 


2.79±0 


16 


-0.87±0.59 


32.6±16.8 


104.3±7.3 


unclassifiablc 


26 


13.2 


12 


3 


31.07±1.40 


7.27 





71 


0.52 


0.043 


1 


521 


2.07±0 


18 


-1.12±0.53 


78.3±16.9 


105.1±8.0 


(0.48) A/BS? 


27 


13.8 


18 


9 


30.00±0.91 


7.67 





80 


0.35 


0.038 


1 


515 


1.36±0 


15 


-1.81±0.43 


90-3±08.9 


78.5±6.5 


(0.94) BHB 


28 


15.2 


16 





29.61±1.09 


7.35 





73 


0.41 


0.043 


1 


512 


2.68±0 


16 


-0.74±0.49 


88. 8± 11.5 


64.9±4.0 


(0.68) BHB 


29 


15.9 


16 




34.59±1.06 


8.57 





77 


0.44 


0.033 


1 


523 


1.83±0 


16 


-0.69±0.54 


92.3±12.7 


76.8±4.8 


(0.51) BHB 


30 


10.5 


11 


2 


25.17±1.34 


6.22 





76 


0.53 


0.059 


1 


509 


3.66±0 


18 


-0.04±0.54 


47.4±21.3 


77.0±4.3 


unclassifiablc 


31 


20.1 


19 





43.93±0.97 


11.20 





80 


0.42 


0.030 


1 


530 


1.28±0 


16 


-1.51±0.44 


41.2±11.3 


47.0±6.7 


(0.00) A/BS 


33 


18.2 


16 


3 


38.49±1.09 


10.13 





83 


0.44 


0.038 


1 


525 


2.50±0 


16 


-0.96±0.45 


133.7±12.9 


40.1±5.0 


(0.08) A/BS 


34 


15.0 


20 





31.80±0.77 


8.66 





88 


0.30 


0.036 


1 


516 


0.96±0 


15 


-2.05±0.46 


119.2±12.3 


80.0±6.6 


(0.98) BHB 


35 


11.6 


13 




27.81±1.16 


6.04 





63 


0.44 


0.037 


1 


513 


1.28±0 


17 


-1.66±0.44 


104.5±12.6 


84.5±6.6 


unclassifiablc 



Table 3. Spectroscopic data for the horizontal branch star candidates. 
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(B-V) 

Figure 4. Call K line (3933 A) EW(A) for the sample plotted 
against (B — V)o. The curves represent lines of constant metallic- 
ity for [Fe/H] = -1.0, -2.0 and -3.0 taken from Wilhelm et al. 
(1999). The straight line represents a best fit to stars in the 
Pleiades and Coma clusters assumed to be of solar metallicity. 
The vertical line at (B — V)o = 0.05 is the limit for which metal- 
licities can be determined. The 8 stars classified as BHB stars are 
marked by filled circles, and the 12 stars classified blue straggler 
are marked by open circles. 



In Figure 03 we plot the two diagnostic diagrams for 
the 20 classifiable stars in the survey. The two figures are 
explained as follows. Figure EI a ) shows the Do. 15 — colour 
method. The average values of Do. 15 for H7 and H<5 against 
(B — V)o are plotted for the 20 candidates. In Paper I we 
showed that reliable classification by this method requires 
the uncertainty on (B — V)o to be less than 0.03 mag. Un- 
fortunately this is untrue for most of the stars in our sample 
(Table 1). For this reason we give this method lower weight 



in the final classification. Figure^Jb) shows the Scale width- 
Shape method. The line-profile quantities b and c, averaged 
for H7 and H5 are plotted. The solid lines show the classifi- 
cation boundaries, from Paper II, with high-surface gravity 
stars (i.e. main-sequence A stars or blue stragglers, here- 
after A/BS) above the line, and low-surface gravity stars 
(i.e. BHB stars) below the line. In both plots stars classified 
BHB are plotted as solid symbols and stars classified A/BS 
are plotted open. As we discuss below, the three triangles 
are stars that have ambiguous classifications, i.e stars that 
are classified as BHB by one classification method and not 
the other. 

Inspection of Figure |S] provides the following informa- 
tion. Of the 20 candidates, eight are classified BHB by 
the Do. 15 — colour method. The Scale width-Shape method 
classifies seven stars as BHB. There are six stars classified 
BHB by both methods. A total of nine stars are classi- 
fied BHB by one or other of the methods. There is clearly 
close agreement between the two classification methods, but 
there are three stars with ambiguous classifications. Before 
considering these further, we note that the three radial ve- 
locity standards (Table 2), previously classified BHB from 
high-resolution spectroscopy, are all unambiguously classi- 
fied BHB in both plots. 

The uncertainties on each parameter define the 2D 
probability distribution functions for any point. By integrat- 
ing these functions below the classification boundary we can 
compute a probability P(BHB) that any star is BHB. We 
can then average the probabilities for the two classification 
methods, to improve the classification. We have computed 
these probabilities for each star, giving twice the weight to 
the Scale width-Shape method when averaging (because, as 
mentioned above, the Do. 15 — colour method is affected by 
the relatively large colour errors). As in previous papers, 
stars with P(BHB) > 0.5 are then classified BHB, and 
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Figure 6. Classification probability against V ga i for all the BHB 
candidates. Large circles represent the 20 stars with spectra of 
high S/N. The BHB stars which have P > 0.5 are shown by 
filled symbols and A/BS stars are marked by open symbols. Also 
plotted, as small symbols, are the nine stars with unreliable classi- 
fications. Finally the five unclassifiable stars are marked by small 
crosses, at P = 0.5. 
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-62.3 


65.1 


27 


310.805 


63.473 


20.8 


76.6 


28 


312.874 


63.489 


21.9 


62.9 


29 


313.702 


62.623 


24.2 


74.7 


34 


315.536 


61.130 


49.9 


77.6 


f225-0011 


290.267 


61.822 


-27.0 


65.0 


f241+0237 


298.280 


65.159 


-32.0 


66.0 


1249+0146 


303.159 


64.372 


-133.0 


53.0 


05 


262.690 


56.439 


50.7 


99.4 


11 


289.614 


63.028 


-112.7 


101.4 



Table 4. Summary of kinematic information for the BHB and 
carbon stars that are associated with the stream. The two BHB 
stars below the line do not belong to the stream. 



stars with P(BHB) ^ 0.5 are classified A/BS. Based on 
the Monte Carlo simulations of Paper I, we would expect 
the sample of BHB stars defined in this way to be contami- 
nated by A/BS stars at no more than the 10% level, which 
we consider satisfactory. Col. (14) of Table 3 provides the 
averaged probabilities, and corresponding classifications, for 
the 20 classifiable stars. 

Of the three stars with ambiguous classifications, 2 are 
classified BHB. The third star classified A/BS, star 9, is the 
object with the smallest value of D0.15, and the smallest 
value of c. The small value of c indicates a colour substan- 
tially redder than (although compatible with) the measured 
value of (B - V) = 0.16 ± 0.04. 

We also provide the classification probabilities for the 
nine stars with inadequate spectroscopic S/N, but for these 
the classifications are given as BHB? or A/BS? to indicate 
that they are not reliable. Finally the five stars with EW 
H 7 < 13A are labelled unclassifiable. 

4.2 Velocity dispersion 

Table 2] contains a summary of the kinematic properties of 
the final sample of eight BHB stars. Listed there are the 
Galactic coordinates / and b, and the Galactocentric radial 
velocity and distance, V 9a ; and r respectively. To convert 
the heliocentric quantities to Galactocentric quantities, the 
heliocentric radial velocities are first corrected for solar mo- 
tion by assuming a solar peculiar velocity of (U, V, W) = 
(-9,12,7), where U is directed outward from the Galactic 
Centre, V is positive in the direction of Galactic rotation at 
the position of the Sun, and W is positive toward the North 
Galactic Pole. We have assumed a circular speed of 220 km 
s _1 at the Galactocentric radius of the Sun (Rq — 8.0 kpc). 
Table 0] then, distills the main observational result of the 
paper, a sample of distant BHB stars with measured radial 
velocities. The Table also includes three carbon stars, desig- 
nated by their coordinates, which are introduced in Section 
5. 

We find, after quadratically subtracting the measure- 
ment errors in the same manner as Norris & Hawkins (1991), 
that the measured dispersion of the radial component of the 
Galactocentric velocity dispersion for our BHB sample is 



58±15km s _1 at a mean heliocentric distance of 80 kpc. In 
Table |S] we compare this value against the measured veloc- 
ity dispersion of a variety of samples. The sample of remote 
BHB stars is referred to as Sample A, and listed in the first 
line of Tabled] The first comparison sample, Sample B, com- 
prises the 60 BHB stars 11 < R < 52 kpc, mean R — 28 kpc, 
from Paper II, which is the largest sample at such distances. 
Sirko et al. (2004a) have also isolated large samples of dis- 
tant BHB stars using the SDSS. They split their sample into 
a bright (g < 18) subsample, which is contaminated by blue 
stragglers at the level of about 10% (i.e. similar to the work 
presented here), and a faint subsample (g > 18), which is 
contaminated at about 25%. If we consider only their clean 
bright sample, here Sample C, then a — 99.4 ± 4.3kms~ 
(Sirko et al. 2004b), at mean distance 16 kpc. Sample D con- 
sists of the 12 stars in Table 3 classified BS, with measured 
velocity dispersion 129 ± 26 km s -1 , at mean distance 40 kpc. 
Finally considering the Galactic satellites discussed in §1, se- 
lecting the nine satellites within the distance range of our 
remote BHB sample, i.e. 65 < R < 102 kpc, Sample E, we 
measure a velocity dispersion 134±32 km s _1 , at a mean dis- 
tance 82 kpc. The velocity dispersions, mean distances, and 
sample sizes, of these four samples are entered in columns 
2 — 4 of Table |^| In the final column we list the probability 
that the measured value for the remote BHB stars could be 
drawn from the same population as each of the four compar- 
ison samples, as measured by the F-test. At better than the 
95% significance level, our sample of remote BHB stars has 
smaller velocity dispersion than BHB stars at much smaller 
radii (16 kpc, Sample C), BS stars of similar apparent mag- 
nitude at intermediate radii (40 kpc, Sample D), and satel- 
lites at the same radii (82 kpc, Sample E). Compared to the 
sample of BHB stars at 28 kpc, the difference is not signifi- 
cant. We conclude that the velocity dispersion of the remote 
BHB stars is anomalously low, and in the following section 
we seek an explanation. 



5 DISCUSSION 

In considering the anomalously low velocity dispersion of 
the remote BHB stars, we first check that the significantly 
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Sample 


Population 


Dispersion 


< R > 


Size 


F-test 






[km s _1 ] 


[kpc] 




[prob.] 


A 


BHB 


58 ± 15 


80 


8 




B 


BHB 


108 ± 10 


28 


(50 


0.09 


C 


BHB 


99.4 ±4.3 


16 


733 


0.02 


D 


BS 


129 ± 26 


40 


12 


0.02 


E 


satellites 


134 ± 32 


82 


9 


0.04 



Table 5. Comparison of measured velocity dispersion of a variety of populations in the Galactic halo, against the measured velocity 
dispersion of our new sample of eight remote BHB stars at mean distance (< R >) of 80 kpc, Sample A. 



different velocity dispersion between the BHB and BS stars 
is robust. In Fig.||j]we plot classification probability against 
Vgai- In this plot, the large symbols represent the 20 stars 
with spectra of high S/N. BHB stars have P > 0.5, and, as 
usual, are shown by filled symbols, with A/BS stars marked 
by open symbols. Also plotted, as small symbols, are the nine 
stars with unreliable classifications. Finally the five unclas- 
sifiable stars are marked by small crosses, at P = 0.5. This 
plot shows a clear difference in the kinematics of stars at the 
bottom of the diagram (high probability A/BS, large veloc- 
ity spread), and at the top of the diagram (high probability 
BHB, small velocity spread). There is no evidence for any 
BHB stars with large values of |V ga i| that have been missed, 
because they fall just outside the classification boundary, or 
because they were unclassifiable because the spectra are of 
low S/N. Therefore, the difference in velocity dispersion be- 
tween the two populations is quite robust to the method of 
classification. 

Another concern we had was the possibility that the 
sample is contaminated by misclassified blue stragglers in 
the Sagittarius stream. Fig. 2 plots a against g* of the ini- 
tial list of candidate BHB stars. As we discussed earlier this 
colour-selected sample of candidates is expected to be con- 
taminated by blue stragglers, because of the large photo- 
metric errors at these faint magnitudes, and this is appar- 
ently confirmed by the high-density of stars at a > 200°, 
where presumably most of the stars are blue stragglers. It 
was therefore worrying that the eight BHB stars, marked in 
the upper diagram by filled circles, mostly lie close to the 
boundary a = 200° in this plot, whereas one might expect 
them to be more uniformly scattered over the RA range. 
If these stars are misclassified blue stragglers, on the edge 
of the Sagittarius stream, this would provide a natural ex- 
planation for the small velocity dispersion. However, if there 
are any Sagittarius blue stragglers a < 200° in our candidate 
list, most will be classified blue straggler, and the reduction 
in velocity dispersion would be greatest in our A/BS sample 
- the opposite of what is seen. 

We conclude from the foregoing discussion that we have 
succeeded in defining samples of BHB and A/BS stars, with 
small contamination, that show distinct kinematic proper- 
ties. Indeed the fact that the measured velocity dispersions 
of the two populations are significantly different is confirma- 
tion of the reliability of the classification methods. In seeking 
an explanation for this difference, a number of possible dy- 
namical explanations could be pursued, for example that the 
stellar orbits change from being predominantly radial to pre- 
dominantly circular at large radii (e.g. Sommer-Larsen et al. 
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Figure 7. Upper: Plot of V ga i versus RA for the six BHB stars 
(shown as solid squares) and three carbon stars (solid triangles) 
that are potential members of a stream. The plot suggests that 
these stars are located near the turning point of an orbit (where 
the radial velocity changes sign). The turning point is located at 
an RA of about 195°. Lower: A plot of RA versus Dec. for the 
stars in the stream. The filled symbols are the stars with V ga ; > 
and the open symbols are for V ga i < otherwise the symbols 
are the same as in the upper plot. The BHB stars are confined to 
the northern stripe of the EDR, data set shown by the horizontal 
lines. 



1997). While this is possible, a more convincing explanation 
was immediately apparent. 

Six of the BHB stars are confined to a small region 
of space with 190° < a < 200°, 63 < r < 78 kpc (and a 
small range in <5). The average distance of these stars from 
the centre of the Galaxy is 70.6 kpc. These stars are there- 
fore confined to a very small fraction of the volume sur- 
veyed. The velocity dispersion of these six stars (corrected 
for measurement errors), 42 ± 12km s _1 , is too large to asso- 
ciate them with a bound object i.e. a low surface-brightness 
dwarf galaxy. The key to understanding the anomalous ve- 
locity dispersion is provided in Figure [Tfupper), which plots 
Galactocentric radial velocity Vg a i versus RA for these six 
BHB stars, marked by squares. A correlation is evident, in- 
dicative of streaming motion, perhaps associated with a dis- 
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rupted satellite. Figure |T|lowerl plots the position on the 
sky of the six stars. To investigate further the possibility of 
a stream, we searched the catalogue of faint carbon stars of 
Totten & Irwin (1998) in the vicinity, confining ourselves to 
the ranges RA 160° < a < 200°, Dec. ± 5°, and to similar 
distances (improved distance estimates taken from Totten, 
Irwin & Whitelock, 2000). Three stars meet these criteria, 
and have been added to both diagrams in Figure |7| marked 
as triangles. Remarkably the three stars appear to add to 
the evidence of a stream. Details of the three stars are pro- 
vided in Table [I] (The velocity errors, taken from Totten, 
Irwin & Whitelock (2000), are 4, 4 and 6 km s _1 in the order 
the stars appear in the Table) 

The correlation between V^ a i and RA evident in Fig- 
urcEIupper), encompasses V ga i = 0, which would correspond 
to a turning point in the orbit at RA of ~ 195° . In order to 
investigate this trend in more detail, we consider orbits in 
the spherical potential 



*(r) 



GM 



Ion 



V a 2 + r 2 + a 



(8) 



The scale length a = 178.0 kpc and the mass M = 
2.0 x 10 12 M are chosen to match those estimated for the 
halo of the Milky Way (Wilkinson et al. 2003). We inves- 
tigate whether the trend in Figure [7| can be reproduced by 
a plausible Galactic orbit as follows. First, we assume that 
the orbit has a turning point in the RA range 190 — 200° 
and in the distance range 50 — 80 kpc. We then choose val- 
ues for the line-of-sight distance do, Galactic latitude bo and 
longitude lo of the turning point and the values of the two 
components of velocity transverse to the line of sight, Vb,o 
and vi t o. From each set of initial conditions, we integrate the 
orbit in the RA range 190 - 200° and determine V^\{RA). 
Assuming Gaussian errors at on the individual radial veloc- 
ities (and neglecting any errors in the RA measurements), 
the probability that the data (V ga \,i, RAi) were drawn from 
the relation V sa i(RA) is given by 



P{vi os ,i, RAi\l , b , r , «b,o, «l,o) = 



n 



: exp 



(vios(RAi) - uios.i) 2 



2a? 



(9) 



We use a downhill simplex algorithm (the routine amoeba 
in Press et al. 1992) to maximise this probability over the 
five dimensional parameter space. The resulting V\ os (RA) 
relation is shown in Figure|7| The orbit we obtain is strongly 
radial and has an apocentre of 67.1 kpc and pericentre of 
7.5 kpc. The energy and angular momentum of the orbit are 
E = -8.1 x 10 4 (kms _1 ) 2 and L 2 = 1.2 x 10 7 (kpc kms" 1 ) 2 . 
The velocity dispersion of all nine stars relative to this orbit 
is 15 ± 4 km s~ . 

Up to this point, we have not made use of the estimated 
distances to our tracer stars, due to their relatively large un- 
certainties. We can include them in the orbit determination 
in a straightforward manner by multiplying the probability 
in eq @ by a second Gaussian P(d\ os ) given by 

P(di os ,i, RAi\l , b , r , u b ,o, «i,o) = 



: exp 



(d loa (RAi) - d u 



2*L 



(10) 



Here, <Jd,i is the uncertainty in the line of sight distance 



to the ith star. In fact, the inclusion of this term results in an 
almost identical orbit, the distance uncertainties rendering 
the distance estimates of little value in the determination of 
orbital parameters. 

We initially analysed the data using the technique of 
Lynden-Bell & Lynden-Bell (1995). These authors note that 
if a group of stars lie on the same orbit in an assumed spher- 
ical potential, they share the same energy E and total an- 
gular momentum L 2 . Since the energy is given by 



1 



1L 2 



\E\ = *(r) - -vi — , 

II v > 2 r 2 r 2 



(11) 



where u r is the Galactocentric radial velocity, L is the mag- 
nitude of the total angular momentum vector and \P(r) is 
the potential, the stars in a stream lie on a straight line in 
a plot of |-E"r| versus r~ where 



= *(r) - -v 2 

\ I 2 r 



(12) 



If we apply this technique to our data, we obtain or- 
bits with very significantly higher orbital angular momenta 
than that for the orbit which fits the V ga ,i versus RA relation 
above. After consideration of the propagation of the obser- 
vational errors into the (|i?r|, l/?" 2 ) plane, we noted that the 
orbital angular momentum was in fact being governed by 
the orientation of the extended error ellipses caused by the 
large distance uncertainties. In addition, the derived orbit 
does not reproduce the trend of V ga i with RA seen in Fig- 
ure We conclude therefore that this technique can only 
yield useful information about an orbit if the magnitudes of 
the distance errors are significantly smaller than the radial 
range covered by the survey - this is not the case for our 
present sample. 

Given the proximity of the Sagittarius (Sgr) stream, 
it is reasonable to suppose the BHB stars might be part 
of a more distant passage of the stream around the Milky 
Way (e.g. Helmi & White 2001; Dohm-Palmer et al. 2001). 
Comparing the positions and heliocentric velocities of our 
sample of stars to the simulations plotted in Dohm-Palmer 
et al. (2001; their figure 3) reveal that our sample largely 
lies either between or beyond two wraps of the Sgr stream. 
In fact, two of our six BHB stream stars (numbers: 34 and 
28) are inconsistent with the model simulations. More re- 
cently, Law, Johnston & Majewski (2005) produced models 
of the Sgr tidal tails using test particle orbits and iV-body 
simulations in a variety of potentials. Before we are able to 
compare our data with these models we need to convert our 
coordinates to the system defined in Majewski et al. (2003). 
In this coordinate system the zero plane of the latitude co- 
ordinate Bq coincides with the best-fit great circle defined 
by the Sgr debris, as seen from the Sun; the longitudinal 
coordinate Aq is zero in the direction of the Sgr core and 
increases along the Sgr trailing stream. Our sample resides 
in the region 249° < A Q < 275° and 7° < B < 20° . At 
these coordinates the whole sample of stars in Table 0] re- 
sides at larger distances than the models predict. Clearly, 
the observation of a larger sample of remote BHB stars in 
the Galaxy halo, along different lines of sight, is essential 
to confirm the reality of the stream. Additionally, we need 
to establish whether the small velocity dispersion measured 
for the eight distant BHB stars discovered in this paper is 
actually because six of the stars are associated with a co- 
herent structure, or because the velocity dispersion of the 
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whole population of outer halo BHB stars falls steeply with 
radius. However, we note that the apparent dominance of 
streaming motion in our BHB sample lends support to the 
claim of Majewski (2004) that the non-uniform kinematics 
of outer halo K-giants are consistent with that population 
having derived almost completely from accretion. 

We close with a summary of the main points of this 
paper. We have presented the results of a survey of remote 
halo A-type stars selected from the SDSS. Spectroscopy of 
the A-type stars obtained with the VLT produced a sam- 
ple of 20 stars with data of suitable quality for classifica- 
tion into the classes BHB and A/BS. The final sample (Ta- 
ble 4) comprises eight stars classified BHB, at distances of 
65 — 102 kpc from the Sun (mean distance 80 kpc), with he- 
liocentric radial velocities accurate to 12 km s _1 , on aver- 
age, and distance errors < 10%. This is the most distant 
sample of Galactic stars with measured radial velocities, of 
this size. Of the eight remote BHB stars, we find that six 
show a strong trend in V ga .\ with RA, and are consistent 
with a single orbit in a spherical halo potential. The mea- 
sured dispersion of the radial component of the Galactocen- 
tric velocity for this sample is 42 ± 12km s -1 . This value 
is significantly smaller than values measured for samples of 
stars at smaller radii, and for satellites at similar radii. This 
evidence is supported by the existence of three previously 
identified carbon stars with the same kinematics. A sim- 
ple model shows all the stars lying on an orbit with energy 
and angular momentum of E = —8.1 x 10 4 (kms -1 ) 2 and 
L 2 — 1.2 x 10 7 (kpc kms -1 ) 2 . The velocity dispersion of the 
nine stars is 56±13km s~ ; the dispersion relative to the 
calculated orbit is 15±4km s -1 . We conclude that we find a 
strong indication of the presence of a stream but further ob- 
servations are required to trace the full extent of this stream 
on the sky. 
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